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Abstract: An extension of the adjoint SU(5) model with a flavour symmetry based 
on the Z4 group is investigated. The Z4 symmetry is introduced with the aim of leading 
the up- and down-quark mass matrices to the Nearest-Neighbour-Interaction form. As a 
consequence of the discrete symmetry embedded in the SU(5) gauge group, the charged 
lepton mass matrix also gets the same form. Within this model, light neutrinos get their 
masses through type-I, type-Ill and one-loop radiative seesaw mechanisms, implemented, 
respectively, via a singlet, a triplet and an octet from the adjoint fermionic 24 fields. It 
is demonstrated that the neutrino phenomenology forces the introduction of at least three 
24 fermionic multiplets. The symmetry SU(5) x Z4 allows only two viable zero textures 
for the effective neutrino mass matrix. It is showed that one texture is only compatible 
with normal hierarchy and the other with inverted hierarchy in the light neutrino mass 
spectrum. Finally, it is also demonstrated that Z4 freezes out the possibility of proton 
decay through exchange of colour Higgs triplets at tree-level. 

Keywords: Discrete and Finite Symmetries, GUT, Neutrino Physics 



ArXiv ePrint: arXiv:1303.XXXX 



Contents 

1 Introduction 1 

2 The model 3 

3 Unification and Proton Stability 7 

4 Effective Neutrino Textures 11 

5 Numerical Results 14 

6 Conclusions 17 
A Matter and Higgs representations 19 
B The Potential 20 
References 21 



1 Introduction 

Grand Unified Theories (GUT) are natural extensions of the Standard Model (SM) and 
provide an appealing framework for the search of the theory of flavour. Most GUT models 
try to unify the three gauge couplings of SM in a unique coupling within a simple group. 
This is sustained by the fact that the SM gauge couplings seem to unify at high scale, A 
-|^q15-17 QgY^ when they evolve through the renormalisation group equations. In such GUT 
constructions, not only the SM gauge coupling unify, but also the SM fermions are tight 
in larger multiplets opening the possibility for the implementation of a flavour symmetry. 
Another important signature of most GUTs is the prediction for proton decay [1], which 
has not yet been observed and severely constrains these models. 

The first GUT model was realisable within the SU(5) gauge group [2] in 1974. This 
minimal model fits the fifteen SM fermionic degrees of freedom in two unique representa- 
tions: 5* and 10 , per generation. It is well established that this model is ruled out, since 
it does not reproduce the correct mass ratios among the charged leptons and down-type 
quarks and also the particle content does not lead to an accurate gauge coupling unification. 
During the last decades, many attempts have been proposed in the literature in order to 
construct consistent GUT models [3-6] based on the SU(5) group. In particular, the mass 
mismatch between the charged leptons and down- type quarks in the minimal SU(5) can be 
easily corrected if one accepts higher dimension operators in the model without enlarging 
the field content [3, 5]. Alternatively, one can build a non-renormalisable solution where 
the mass mismatch is explained by adding an extra 45 Higgs multiplet [7, 8]. 
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Although GUT multiplets contain both quark and lepton fields this is not enough to 
fully determine the properties of their observed masses and mixings. Indeed, GUT models 
do not solve the "flavour puzzle" present in the SM, however the new GUT relations 
among quark and lepton Yukawa matrices are an excellent starting point for building a 
flavour symmetry. There have been many approaches to understand the intricate "flavour 
puzzle" in the context of GUTs. An attractive possibility is to assume the vanishing of 
some Yukawa interactions by the requirement of a discrete symmetry, so that new "texture 
zeroes" appear in the Yukawa matrices [9-14]. Symmetries may predict new relations 
among fermion masses and their mixings. Nevertheless, the opposite is not true in general, 
because zeroes in the Yukawa matrices can also be obtained by performing some set of 
transformations (weak basis transformations) leaving the gauge sector diagonal [15-17]. 

The Nearest-Neighbour- Interaction (NNI) is an example of a weak basis in which the 
up- and down-quark mass matrices, and M^, share the same texture-zero form: 



where the constants Au^d > d ' ^u,d , B'^ d Cu,d are independent and complex. Since 
this parallel structure is a weak basis, no physical predictions can be made unless further 
assumptions are considered. This is the case of the Fritzsch Ansatz [18, 19] where, in 
addition of the NNI structure, one requires Hermiticity for both and M^, which cannot 
be obtained through a weak basis transformation. It is well known that the Fritzsch 
Ansatz can no longer accommodate the current experimental quark mixings. However, it 
was shown in ref. [20] that deviations from the Hermiticity around 20% were compatible 
with the experimental data. 

It was shown in ref. [20] that it is possible to obtain up- and down-quark mass matrices 
Mu and with the NNI structure through the implementation of an Abelian discrete 
flavour symmetry in the context of the two Higgs doublet model (2HDM). In that context, 
the minimal realisation is the group Z4. In a general 2HDM, a NNI form for each Yukawa 
coupling matrices cannot be a weak basis choice. Indeed, the requirement of the Z4 sym- 
metry does imply restrictions on the scalar couplings to the quarks, although one gets no 
impact on the quark masses and the Cabibbo-Kobayashi-Maskawa matrix [21, 22]. 

The goal of this article is to study whether it is possible to implement a Z4 flavour 
symmetry, as in refs. [20, 23], that leads to quark mass matrices and with the 
NNI form in the context of the adjoint-SU(5) Grand Unification [6]. The requirement of a 
flavour symmetry that enforces a particular pattern in both up- and down-quark Yukawa 
couplings has phenomenological implications for the leptonic sector. The adjoint-SU(5) 
model strengths the weak points left in ref. [23], namely it improves the unification of the 
gauge couplings, solves the mass mismatch between the charged leptons and down-type 
quarks at renormalisable level, alleviates the constraint imposed by the proton lifetime and 
introduces a richer mechanism to generate light neutrino masses. 

The minimal version of the adjoint-SU(5) model [6] consists in adding to the minimal 
SU(5) version a 45 Higgs representation, 45ff , together with an adjoint fermionic field, /j(24) 
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to generate light neutrino masses. In this minimal setup one /o(24) is enough to account 
for the observed low-energy neutrino data compatible with a non-zero lightest neutrino 
mass [24-27]. Nevertheless, the nature of the Z4 symmetry requires an extension of the 
number of p(24) fields, 7224 . The fermionic p{2A) fields make possible the generation of 
neutrino masses through the interaction with the 45// field, since gauge interactions forbid 
a singlet right-handed neutrino to couple to 45//. Neutrino masses arise from three different 
types of seesaw mechanisms: type-I [28-31], type-Ill [32, 33] and radiative seesaw [34-36]. 
The radiative seesaw is realizable through the octet-doublet 'S'(8,2) of the 45// multiplet at 
the one-loop level. 

This paper is organised as follows: in section 2 the SU(5) x Z4 model is described in 
detail. Next, in section 3 we address the issues of the unification of the gauge couplings 
as well as the phenomenology of proton decay in the model. The successful textures for 
the effective neutrino mass matrix together with some comments on the generation of the 
baryon asymmetry of the universe through leptogenesis are discussed in section 4. Then, 
our numerical results showing the viability of the leptonic textures considered are sketched 
in section 5. Finally, the conclusions are drawn in section 6. 

2 The model 

The adjoint-SU(5) model [6] contains three generations of 5* and 10 fermionic multiplets 
which accommodate the SM fermion content. In addition one adjoint fermionic multiplet 
/9(24) is introduced for the purpose of generating the light neutrino masses and mixings. In 
this minimal version, three non-vanishing light neutrino masses arise from three different 
seesaw mechanisms, as it will become clear later. Since our aim is to enlarge the symmetry 
of the Lagrangian with an extra Z4 symmetry it forces us to consider n24 copies of 24 
fermionic field. 

The Higgs sector is composed by an adjoint multiplet, S(24), a quintet 5// and one 
45-dimensional representation 45//. Details of the SM fields contained in the GUT repre- 
sentations are given in appendix A. 

The adjoint field S has the usual role to break spontaneously the GUT gauge group 
down to the SM group, i.e. SU(3)c x SU(2)l x U(1)y, through its vacuum expectation 
value (VEV), 

(S) = ^diag(2,2,2,-3,-3). (2.1) 
VoU 

The Higgs quintet and the Higgs 45-plet give rise to two doublets. Hi € 5// and H2 G 45//, 
at low-energies and two massive SU(3) colour triplets, Ti G 5// and T2 € 45//. It is 
essential that the triplets Ti^2 have masses around the unification scale while the doublets 
Hi^2 should remain at the electroweak scale in order to prevent rapid proton decay - the 
so called doublet-triplet splitting problem. The representations (3,1,4/3) and (3,3,-1/3) 
in 45// could also induce proton decay. However, it has been shown that the former 
representation does not contribute at tree level to proton decay, while some states of the 
latter representation contribute with a constraint milder than the one given by the triplet 
T2 [37]. 
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Many mechanisms were proposed in order to avoid the doublet-triplet-splitting prob- 
lem. One possibility that can be easily invoked within this framework is the missing partner 
mechanism [38, 39], which consists in having the bosonic representations 50, 50* and 75 
instead of the adjoint S to break the GUT group. In the missing partner mechanism the 
scalar doublets are naturally massless. 

The role of the scalar fields 5^ and 45h is then to break the SM group to SU{3)c x 
U{l)e.m. through their VEVs 

{5h)^ = iO,0,0,0,vs) , (2.2) 

and 

(ASnf) = v,s{s^a-'^S2^t) , «,/3 = l,...,4, (2.3) 

that are related as 



|f5|V24|w45p = (\/2Gf) ^ = (246.2 GeV)^ (2.4) 



where Gp is the Fermi constant. 

Since all representations are well defined, we can specify the nontrivial transformations 
of each bosonic/fermionic field R under the Z4 fiavour symmetry as: 

R' = e^'-r R, Q{R)eZ4. (2.5) 

The purpose of the discrete symmetry is to obtain the quark mass matrices, M„, Md, 
with the NNI form at low energy scales. In order to preserve the Z4 symmetry below the 
unification scale it is required that Q(S) = 0. In order to implement the missing partner 
mechanism, the extra bosonic fields, 50, 50* and 75 , must also be trivial under Z4 . Thus, 
below the unification scale A the Z4 group is preserved in higher orders of perturbation 
theory, provided that no Nambu-Goldstone boson appears at tree-level due to an accidental 
global symmetry [40]. At low energies one obtains a two Higgs doublet model with extra 
fermions invariant under Z4, which gets broken once the doublets acquire VEVs. 

The Z4-charges are assigned as follows. First we make the choice that the Abn couples 
to the bilinear IO3 IO3 , which implies that 

Q(45h) = -2Q(103). (2.6) 

This particular choice does not eliminate any texture on the leptonic sector obtained when 
varying the fermionic Z4-charges. Thus, the most general fermionic Z4-charges that lead 
to NNI for the quark mass matrices Mu^d are 

Q(lOi) = (3g3 + 0, -93 - 93) , 
Q(5^) = (^3 + 20, -3(?3, -93 + </>), 

where (j) = Q{5h) and = Q(103). The charges for the 7^24 adjoint fermions are left free 
and only some combination of them will lead to realistic effective neutrino mass matrices, 
as we will see. 
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In this model, the most general Yukawa interactions are given by the following terms: 



+ \,i Tr {p, pi S) + (r^)^^ 5*, (pkTp {buf + {Tl)^^ 5*^ (p,)} {45H)f + H.c 



(2i 



where a, /3, • • • = 1, . . . , 5 are SU(5) indices, i,j are generation indices and k,l = 1, . . . , 7124 . 
Notice that the Yukawa matrix and A as well as the mass matrix M are symmetric 
while is antisymmetric. Taking into account the charges given in eq. (2.7), the Yukawa 

1 2 

coupling matrices Fj^ are given by. 





(2.9a) 



(2.9b) 



b', 



The up- and down-quark masses as well as the charged lepton masses are given by 

Mu = Av5ri + 8^745 



(2.10) 



Substituting the eqs. (2.9) in the mass matrices given in eqs. (2.10) one concludes that all 
the matrices M^^d^e share the NNI structure. The up-quark mass matrix is no longer 
symmetric and the mismatch between the down-type and charged lepton matrices is now 
explained as: 



(2.11) 



The mass matrices of the fermions poi P3 ^^nd pg arising from the fermionic-24 fields 
are given by: 



Mo 

M3 
Ms 



M - 
M - 
M + 



a 

\/30 
3cr 



^30 
2a 



30 



(2.12) 



Due to the fact that the Higgs field S is trivial under Z4, the matrices M and A share 
the same form and this is also valid for the Majorana matrices Mo,3,8- From the Yukawa 
interactions written in eq. (2.8), one can infer the Yukawa couplings for the pQ , /J3 and pg 
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Figure 1. Different seesaw mechanisms present in the model. 



fermion fields, which are then given by 



Y 



COS a 



+ 



2V2 L 5 



2V2 



sma 



kl 



ljia2 poi H' 



+ [cosa (r^)^; -3sina(r^)^J lJia2P3iH (2.13) 
--^ [sma{Tl)j^^ + 3cosQ(r^)^J llia2p3iH' 

-^^llia2TT{S^S,2)P8l) , 



where »S'(8,2) is the scalar octet-doublet belonging to the 45^ representation and the doublet 
space {Hi,H2)^ has also been rotated in terms of new doublets {H, H')^ such that {H) = v 
and {H') = by the appropriate transformation: 



H \ I cos a sin a\ I Hi 
H' ) I — sin a cos a ) \H2 



(2.14) 



with tan a = V45 / V5 . 

Taking into account the Majorana mass matrices Mo,3,8 and the Yukawa interactions 
given by eq. (2.13) one can derive, through the seesaw mechanism, the effective neutrino 
mass matrix mj, which receives three different contributions, as drawn in figure 1. One has 
type-I seesaw [28-31] mediated by the fermionic singlets pofc and type-III seesaw [32, 33] via 
the exchange of the SU(2)-triplets ps^. There is still the possibility of generating neutrino 
masses through the radiative seesaw [34, 35], that involves at 1-loop the fermionic psk and 
the scalar doublet-octet ^(g 2) present in the 45//. The neutrino mass matrix obtained after 
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integrating out the fields responsible for the seesaw mechanism reads as 



k=l 

where uiq , are given by 



(2.15) 



mo =^(— r, + smar,j , 



V2 



7715' = (cos a rj, — 3 sin a F^) , 



(2.16) 



and Mg ^ are simply the mass-eigenvalues of the Majorana matrix Mg . The unitary matrix 
Us does the rotation of the yukawa matrix to the basis where the matrix Mg is diagonal. 
The coefficient is a linear combination of the coefficients in the Higgs potential terms 
given in eqs. (B.2c) and (B.2f) in the appendix (B). The loop function F{x) is given by 



If one assumes Mg k ^ Mq ^ > -^3 k then it suppresses the 1-loop radiative seesaw contri- 
bution. 

Before closing this section, it is important to comment about the Higgs potential. The 
most general Higgs potential is given explicitly in the appendix (B). Notice that, terms 
involving simultaneously the fields bu, 24^^ and Abn are forbidden by the Z4 symmetry. 
This gives rise to an accidental global continuous symmetry which, upon spontaneous 
electroweak symmetry breaking, would lead to a massless Nambu-Goldstone boson at tree 
level [41]. A simple way to cure this problem is by adding a complex SU(5) Higgs singlet 
5 non-trivially charged under Z4 , i.e. Q{S) = — 2 (73 — , where its potential is given by 

Vs = (Asb 5; 2A} 45f S + H.c.) + ^ 4 \S\^ + A5 + A'^ (5^ + H.c.) , (2.18) 

and leads at low-energy to an effective interaction, once the scalar S acquires vacuum 
expectation value, 

AsB <T (5) i^J i72 + H.c. , (2.19) 
which softly breaks the symmetry Z4 . 



3 Unification and Proton Stability 

According to the previous section, between the unification scale and Mz = 91.1876 it 
0.0021 GeV scale [38] one has a 2HDM with extra fermions, namely pofc ) PSfc > and pg^ , and 
the two scalars S3 and Eg that can have lower masses. The coulored triplets, T\ and T2 , 
and the other scalars contained in the Abn are set their masses arround the GUT scale. We 
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also assume Msg ~ -^Ss- The running of the three gauge coupling constants (i = 1, 2, 3) 
in the 2HDM with extra particle content can be obtained easily at the one-loop level as 



bl 


log 




2tt 




b2 


log 




2tt 






b3 


log 




27r 





«3-'(,.) - <m) - i ^4 j - E i 10. (^) , (3,10 

where ai = 5/3 

Q^i/i ct2 — and 03 — OLg] the hi constants are the usual one-loop beta 
coefficients corresponding to the 2HDM, listed in section A. M/ denotes an intermediate 
energy scale for extra particle / between the electroweak scale Mz and the GUT scale A, 
and the coefficients hi account for the new contribution to the one-loop beta functions hi 
above the threshold M/. At the unification scale A, the gauge couplings Oj obey to the 
relation 

ajj = ai(A) = Q2(A) = a3(A) . (3.2) 

To get some insight into the unification in the one-loop approximation, let us define 
the effective beta coefficients Bi [42], 

B, = h, + Y,blru (3.3) 

where the ratios < r/ < 1 that takes into account the intermediate scales are given by 

^ ln(A/M,) 

' \n{K/Mz) ^ ' 

It is also convenient to introduce the differences Bij = Bi — Bj, define as 

B.,=B!j^^ + Y.^h-^^ (3-5) 
where Bfj^^^ corresponds to the 2HDM particle contribution and 



Aij^bf-b^. (3.6) 



The following i?-test is then obtained, 



a 



B23 



sin^ dw 

a 



B^-f = ^ ^' (3.7) 



5 5 



together with the GUT scale relation 
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Notice that the right-hand sides of eqs. (3.7) and (3.8) depend only on low-energy 
electroweak data, using the following experimental values at Mz [38] 

= 127.916 ±0.015, (3.9) 

sin2 6'vy = 0.23116 ±0.00012, (3.10) 

a, = 0.1184 ±0.0007, (3.11) 

the above relations read as 

B = 0.718 ±0.003, 

/ A \ (3.12) 
Bi2 In — = 185.0 ±0.2. 
\MzJ 

The coefficients Bij that appear in the left-hand sides of eqs. (3.7) and (3.8) strongly 
depend on the particle content of the theory. For instance, considering the SM particles 
with uh light Higgs doublets, one has bi = 20/5 + nH/lO, 62 = — 10/3±n///6 and 63 = —7, 
so that these coefficients are given by 

Si2 = 36/5, ^23 = 4. (3.13) 

where tlh = 2 is set and B = 5/9 is then not compatible with the calculated value in 
eq. (3.12) and clearly, the B-test fails badly in the 2HDM case, so that extra particles are 
needed. In table 1 we present the relevant contributions Ajj to the B^j coefficients of our 
setup which include, besides the 2HDM threshold, the fermions poj. , p^j. , and psfc > ^nd the 
two scalars S3 and Sg are considered. For simplicity, we assume the remaining particles at 
the unification scale and therefore they do not contribute to the gauge coupling running. 



Table 1. The relevant A^^ contributions to the Bij coefficients in the SU(5) x Z4 model. 





2HDM 


P3 


P8 


S3 




A12 


36/5 


-4/3 





-2/3 





A23 


4 


4/3 


-2 


2/3 


-1 



Notice that eqs. (3.12) require Ajg < and A23 > 0, it becomes clear from table 1 all 
extra particles considered in the running improve the unification. We have scanned different 
ratios rj and we obtained a large range of solutions that lead to a perfect unification 
within the experimental errors. The fact that more than one adjoint fermionic field pk is 
present, it improved the range of intermediate scales Mj consistent with unification. It is 
now difficult to find strong correlations among the intermediate scales Mj. In addition, 
the possible values for the unification scale A can vary many orders of magnitude. For 
illustration, in figure 2 we have drawn the mass spectrum of the extra particles included 
in the running. All the solution obtained are in agreement with a unified gauge coupling 
gu = \J A-Tx ja^ < 1, where in our numerics we obtained 37 . 

Concerning the proton decay some comments are in order. In this model there are 
mainly two different sources for proton decay, namely via the exchange of the lepto-quark 
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Figure 2. Four illustrative examples showing the mass spectrum of the adjoint fermionic fields, 
S3 and Es for different unification scales A. 



gauge bosons X, Y or via coulored Higgs triplets. Proton decay in both scenarios are 
mediated four fermion interactions (dimension-six operators). 

The gauge bosons X, Y become massive through the Higgs mechanism with a common 
mass, My, 

Mv = -9W. (3.14) 

To suppress the X, Y boson proton decay channels, one has necessarily that My ^ mp 
(the proton mass) , that leads to the estimation of the proton decay width as [8] : 

5 

r-a?.^. (3.15) 

Making use of the most restrictive constraints on the partial proton lifetime r(p —t- 71*^6+) > 
8.2 X 10^'^ years [38], one can derive a rough lower bound for the X, Y mass scale My, 

My > 4.1 X 10^^ GeV, (3.16) 

which corresponds a a^^ w 37. Since we assume for the unification scale A ~ My, the 
constraint given by eq. (3.16) determines the scale where the gauge couplings should unify 
(for a recent review see [1]). 

The proton decay through the exchange of Higgs colour triplets Ti , T2 is very sup- 
pressed, since their suppression is proportional to products of Yukawa couplings, and 
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therefore they are much smaller than the gauge couplings. Indeed the contribution of these 
dimension-six operators vanishes at tree- level when the Z4 symmetry is exact [23]. The 
dimension-six operators contributing to the proton decay via the colour triplet exchange 
are given at tree-level by: 



E 

n=l,2 



m2 



-{QiQi){QkLi) + {u're'j){uldf 



(3.17) 



It is then clear from the pattern of the Yukawa coupling matrices and given in 
eqs. (2.9) that the only possible non-vanishing contribution of the dimension-six operators 
given in eq. (3.17) involve necessarily fermions of the third generation. One concludes that 
at tree-level the proton does not decay through the four-fermion interactions described by 
the operators given in eq. (3.17). 



4 Effective Neutrino Textures 

The flavour symmetry present in our model constrains the charges of the fermion fields 
to be of the form in eq. (2.7). Such charge assignment does not imply any restriction in 
the neutrino sector. Hence the charges of the 24 dimensional fermionic representations 
responsible for the neutrino masses are free. In order to analyse the possible patterns for 
the effective neutrino mass matrices, niu, we have considered all the possible values for the 
Z4 charges of the adjoint fermionic fields. 

Searching for the minimal SU(5) x Z4 model, i.e. the model with the minimal matter 
content, we have started by the possibility of having only one extra 24 fermionic repre- 
sentation, as in the SU(5) adjoint original model [6]. However, given the particularities 
of the Z4 symmetry, the neutrino mixing pattern that emerges from this picture is now 
not consistent with the experimental neutrino data [43]. Adding a second 24 fermionic 
representation does not solve the problem and again the predicted neutrino mixing angles 
are not in agreement with neutrino oscillations. 

The situation changes when we consider three 24 fermionic representations. In this 
case we obtain different possibilities for the light neutrino mass matrix rrii, that coincide 
with the matrices found in ref. [23], where a similar Z4 flavour symmetry was imposed in 
the context of three right-handed neutrinos. From the various textures for the effective 
neutrino mass matrix rrii, found in the scan, only two solutions can account successfully 
for the low-energy neutrino data, namely 



m,. 




and 



m, 



i(12) 




(4.1) 



where the index (12) refers to the fact that texture- ^4(12) is a permutated form of the 
texture-A through the permutation matrix P(i2) , 

/o 



^(1 



(12) 



1 



1 






(4.2) 
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isomorphic to the symmetric group 5*3. Thus, the effective mass matrices and m,y 
are related by: 

=P,i,)m;^pT)- (4.3) 

The compatibihty of the above textures with the experimental neutrino data has been 
analyzed in detail in ref. [23] where it was found that just the two-zero textures in A and 
^4(12) are phenomenologically viable. An important result is that texture A is compatible 
only with normal hierarchy (NH) while texture A12 turns to be compatible only with 
inverted hierarchy (IH) in the light neutrino mass spectrum. 

At this point, we would like to remark that the textures given in eq. (4.1) have also been 
studied in the literature previuosly. According to the standard terminology [44] , our allowed 
matrices A and ^(12) would correspond to the ones labelled as A2 and Di, respectively. 
Texture Di has been shown to be either disallowed [44-46] or very marginally allowed [47]. 
However, all of these previous works assume a diagonal charged lepton mass matrix while 
here we are considering a charged lepton mass matrix with NNI form. Therefore the results 
in the literature do not strictly apply to our case. 

Up to here we have only considered the Z4 charges for the 24 fermionic fields such 
as the matrices M and A are non-singular, i.e., |M| 7^ and |A| 7^ and therefore the 
Majorana matrices M^^, x = 0,3,8 are non-singular as well. This condition is mandatory 
in order to derive the effective neutrino mass matrix using the seesaw formula in eq. (2.15). 
However, given the flavour symmetry present in our model, among all the possible config- 
urations one can have |M^| =0. In these cases is not defined and therefore we can 
not use the standard seesaw formula, but instead we should consider the singular seesaw 
mechanism. The singular seesaw, first suggested in the context of a GUT framework [48], 
has been considered in relation to different anomalies in neutrino physics such as the Simp- 
son neutrino [49-51] or the LSND signal [52-55]. In both cases the modified seesaw scheme 
has been used to obtain a neutrino mass spectrum with a singlet/sterile neutrino in the 
energy range between light neutrino masses below the eV and the heavy neutrinos at the 
seesaw scale. 

Here we have analysed the neutrino mass spectrum that emerges from the singular 
seesaw mechanism with two and three 24 fermionic representations. According to our 
calculations, it is not possible to obtain an effective neutrino mass matrix compatible with 
experimental data in any of the cases, since we always get too few light neutrino states. 

In summary, our analysis shows that the symmetry requires the presence of at least 
three fermionic 24 representations and no singular seesaw, i.e., [M^,.! 7^ 0. Since the param- 
eters of the last matrix depend on the Z4 neutrino charges Q(24j), we can conclude that 
the neutrino phenomenology has an impact on the Z4 neutrino charges. 

We present in table 2 the Z4 charge assignment for the fermionic fields that leads to 
the successful textures for the effective neutrino mass matrix A and ^(12), discussed above. 
The textures for the matrices Mx and rriQ^ are also shown. The charges for the fermionic 
fields 10 and 5* as well as the Higgs fields 5h and 45h follow the relations in eqs. (2.6) 
and (2.7). It is worth pointing out that all Dirac neutrino matrices "2-0*3 obtained through 
the scan have four texture zeroes. 
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Table 2. The Z4 fermionic field charges for phcnomenologically viable effective neutrino textures. 



M, 



0,3,8 



Q(240 



'0,3 



2(5*) 2(10,) Q{5h) 2(45h) 



^0 * ON 

^ ^ ^ 

vO * */ 



^0 *N 
0*0 
Oy 



ON 
* 
vO * Oy 



(1,2,3) 



(0,1,3) 



^0 *N 

* * 
yO * 

^* ON 
* * 
^* * 0/ 

^0 * ON 

* * 
^* * 0/ 

^0 *N 

* * 

*/ 



(3,1,0) (0,2,1) 

(1,3,0) (0,2,3) 

(1,3,2) (2,0,3) 

(3,1,2) (2,0,1) 






(1,2,3) 



(0,1,3) 



(2,0,1) (1,3,0) 



(2,0,3) (3,1,0) 



(0,2,3) (3,1,2) 



(0,2,1) (1,3,2) 



Leptogenesis 

In this section we would like to briefly comment about the possibility of having leptogenesis 
in our model. 

As already discussed in refs. [26, 27], the out-of-equilibrium decays of the fermionic 
fields pq and ps in the 24 fermionic representation may generate an asymmetry in the 
leptonic content of the universe. In the presence of sphaleron processes this leptonic asym- 
metry would be partially converted into a baryon asymmetry, explaining the observed 
matter-antimatter asymmetry of the universe. Depending on the mass hierarchy among po 
and p3, the main contribution to the leptonic asymmtery will be dominated by the decays 
of one of them. In principle, our model has enough freedom to have different mass spectra 
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for the fermionic fields and tlierefore, in contrast to the results shown in refs. [26, 27], in 
our case the leptonic asymmetry may be generated by the decay of po (eo) or (es). In 
both cases, the expression for the generated CP asymmetry would be proportional to: 



Despite the specific flavour structure of the matrices, induced by the Z4 symmetry 
(see eq. (2.16)), we have checked that in principle there are no cancelations in the terms 
above and therefore the leptonic asymmetry generated by po and pa can be different from 
zero. It is clear that more accurate predictions about leptogenesis would require further 
calculations, considering the effect of the washout over the initial leptonic asymmetry as 
well as the dynamical evolution of the asymmetry with the solution of Boltzmann equations. 
For the moment, however, our goal is just to show that the model presented here has enough 
freedom in the choice of masses and couplings so in principle it is possible to accommodate 
the CP asymmetry. Further considerations as, for instance, the constraints on the model 
coming from the requirement of a baryon asymmetry consistent with the observations will 
be discussed elsewhere. 

5 Numerical Results 

In this section we analyse the phenomenological implications of the effective neutrino mass 
matrices and rriu in eq. (4.1). Since the flavour symmetry is valid under perturbative 
corrections until the breaking of the electroweak gauge symmetry, the form of the Yukawa 
matrices F^'^, F^'^ and the Majorana mass matrices Mo,3,8 remains unchanged. Thus, one 
can extract the predictions for Mg and rriu and confront them with the observed neutrino 
data at Mz energy scale. The effective neutrino mass matrices obtained mj^''^'^^' , as already 
mentioned, are the same as those analysed in ref. [23]. However, the new measurements of 
the reactor mixing angle ^13 [56] have changed the theoretical picture of the light neutrino 
mixings since then. Therefore, it is worth to revisit the previous analysis in ref. [23] to take 
into account these new bounds. 

Without loss of generality one can write the charged lepton mass matrix. Me, and the 
effective neutrino mass matrices, m[f^ as: 




(4.4b) 



(4.4a) 



/ Ae 

Me = Ke A'g Be 
\0 Bi Ce 



(5.1a) 




(5.1b) 
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where the permutation g = eoi (12) according to the table 2 and the constants ^e,v) Bg^i^, A'^, B'^, Ce,u, 
are taken real and positive. The diagonal phase matrix can be parameterised as 

Ke = diag{e'^\e'^\l), (5.2) 

and the phase (p in eq. (5.1b) cannot be removed by any field redefinition. 

Although the number of the parameters encoded in the pair Me,mu is 12, as the 
number of independent physical parameters experimentally observed at low energy, the 
zero pattern exhibited in eqs. (5.1) does imply new constraints among the independent 
physical parameters, as it will be shown. The PMNS matrix U is given by 

U = Oj Kl , (5.3) 

where the orthogonal matrix Og is the one that diagonalises M^mI as 

{Ke Oe)^ Me (Ke Oe) = diag(me , , rrir) , (5.4) 

while the unitary matrix Ui, diagonalises as 

{Pg U^)'^ ml {Pg Uu) = diag(mi , m2 , ms) . (5.5) 

The knowledge of the low-energy neutrino mixings appears in the literature in terms 
of the parameters 9i2 , 9i3 , 623 and 5 of the Standard Parametrisation (SP) [38], defined 
in terms of PMNS matrix U invariants as 

sm6'i2 ~ 



V^-\Ue3\^' 

sin 6*13 = \Ue3\, (5.6) 
sing23^ , 

Vl-|^e3P' 

and the phase 5 is given by the Dirac-phase invariant, /, 

I = Im {Uf,3U:3Ue2U;2) = ^ COS ^13 sin 2012 sin 2^23 sin 2^13 sin 5 . (5.7) 

Due to the fact that the PMNS matrix is not rephasing invariant on the right, one defines 
the Majorana-type phases, ffj, free of any kind of parametrisation as [57]: 

^tj^^rgiU^iU:^) . (5.8) 

It has been shown in ref. [57] that the PMNS matrix can be fully reconstructed by six 
independent Majorana-type phases from eq. (5.8) taking into account that [/ is a uni- 
tary matrix. The Dirac-type phase 5 can therefore be expressed as the difference of two 
Majorana-type phases: 

I = \Uf,3\\Uf,2\ \Ue3\ \Ue2\ sin ((^^3 - ^^3) . (5.9) 

In what follows we will use the three Majorana-type phases 1P23 , ^P23 ^^'^ ■ 
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Table 3. Neutrino oscillation parameter summary from ref. [43]. For Amfi, sin^ 6*23, sin^ ^13, and 
S the upper (lower) row corresponds to normal (inverted) neutrino mass hierarchy. 



parameter 


best fit 


3cr range 


Am^i [10-^ eV] 


7.62 


7.12-8.20 


|Am2j [10-3 eV] 


2.55 
2.43 


2.31 - 2.74 
2.21 - 2.64 


sm^ 612 


0.320 


0.27 - 0.37 


sin^ 023 


0.613 (0.427) 


0.36 - 0.68 


0.600 


0.37-0.67 


sin^ 013 


0.0246 
0.0250 


0.017-0.033 


5 


O.SOtt 
-0.037r 


0- 27r 



In our analysis we have calculated Og numerically using the charged lepton masses 
given at Mz scale in the MS scheme at 1-loop [58, 59] as 

nie = 0.486661305 ± 0.000000056 MeV , (5.10a) 
= 102.728989 ± 0.000013 MeV , (5.10b) 
m-r = 1746.28 ± 0.16 MeV . (5.10c) 

Concerning the neutrino sector we have used the most recent three neutrino data from the 
global fit of neutrino oscillations in ref. [43]. The best fit values and 3a ranges for the 
neutrino parameters are presented in table 3. 

As in ref. [23], here we have varied all the experimental charged lepton masses and 
neutrino mass differences within their allowed range given in eq. (5.10) and Table 3, re- 
spectively. The mass of the lightest neutrino (mi in NH or in IH) was scanned for 
different magnitudes below 1 eV. To reconstruct the PMNS matrix, we have also scanned 
the free parameters Ae , B^, Dy and the phases ki, K2, defined in eq. (5.1). All the 
remaining parameters are calculated in terms of the former ones. The restriction in this 
scan was to accept only the input values which correspond to a reconstructed PMNS ma- 
trix U that naturally leads to the mixing angles 612 , ^23 and ^13 within their experimental 
bounds presented in Table 3. 

Prom our scan we have found that the mass matrix A in eq. (4.1) is compatible with 
a neutrino mass spectrum with normal hierarchy while the texture ^(12) is compatible 
with inverted hierarchy. The allowed ranges for the lightest neutrino masses are mi = 
[0.353,20.884] x lO'^ eV for NH and = [2.575, 15.335] x lO'^ eV for IH. The presence 
of a massless neutrino as well as a quasi-degenerate neutrino mass spectrum are excluded 
in both cases. For the texture A we have found no significant correlations between sin^ ^13 
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sin\, 0,7 




Figure 3. Plot of sin^ ^13 as a function of ma (left panel) and sin^ 623 as a function of sin^ 612 
(right) in the case of texture mi/'^^' and inverted hierarchy. 

and sin^ ^23 as a function of mi, while in the case of texture ^(12) some correlations are 
found, as shown in figure 3. In fact, this correlation is behind the narrower allowed 
range for IH in comparison with the allowed mi range for NH. We have also verified that 
textures A and ^(12) are not compatible with inverted and normal hierarchies, respectively, 
even when the new limits on sin^ ^13 are considered. 

In figure 4 we plot the effective Majorana neutrino mass characterizing the neutrinoless 
double beta decay amplitude ruee with respect to the lightest neutrino mass, mi in the 
case of and normal hierarchy or ms in the case of and inverted hierarchy. The 

shadowed bands correspond to the generic predictions for niee according to the experimental 
neutrino data at 3a, without any further assumption concerning the origin of neutrino 
masses. If we now restrict our calculations to the adjoint SU(5) x Z4 model presented in this 
article, the allowed regions are reduced to the darker pointed regions. The horizontal lines 
in figure 4 correspond to the mge sensitivity that will be reached by the next generation of 
neutrinoless double beta decay experiments (see for instance ref. [60]). There we see that, 
even if a part of the inverse hierarchy band will be experimentally covered in the next years, 
a sensitivity of around 10-30 meV will be needed in order to probe the effective Majorana 
mass predicted by our model. Accessing to the predicted region for normal hierarchy will 
be even tougher, since a sensitivity of the order of 1 meV would be required, far away from 
the more optimistic scenarios. 

6 Conclusions 

In this work we have studied the adjoint-SU(5) x Z4 model. The flavour symmetry imposed 
in the Lagrangian has the purpose to force the quark mass matrices M„, to have the NNI 
form after the spontaneous electroweak symmetry breaking [20, 23]. Due to the fermion 
content of the adjoint-SU(5), the charged lepton mass matrix lias automatically NNI 
form. In this model the light neutrinos get their masses through type-I, type-Ill and one- 
loop radiative seesaw mechanisms, implemented, respectively, via a singlet, a triplet and 
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0.0001 0.001 0.01 0.1 



m [eV] 

Figure 4. Effective Majorana neutrino mass m^f. as a function of the lightest neutrino mass m 
for normal and inverted neutrino hierarchy, as indicated. The upper band shows the experimental 
sensitivity to be achieved in the next years. 

an octet from the adjoint fermionic fields. The SU(5) x Z4 symmetry does not impose any 
constraint on the adjoint fermionic fields. 

We have shown that the model proposed is in agreement with the current experimental 
data. Neutrino mixings and mass splittings as well as the masses of the charged leptons have 
been used to constrain the possible textures of the effective light neutrino mass matrix. We 
have demonstrated that at least three copies of the 24 are needed in order to fully implement 
the Z4 flavour symmetry and simultaneously account for the experimental neutrino data. 
As shown in table 2 only two zero-textures persist: A and ^{12), which are compatible with 
normal and inverted hierarchies, respectively. 

One of the main phenomenological implications of the model studied is the prediction 
of a hierarchical neutrino mass spectrum not compatible with a massless neutrino. This 
result is particularly important since the neutrino mass spectrum predicted can be used to 
prove or disprove the model in the near future. At present our results are in agreement 
with the constraints coming from neutrinoless double beta decay [60] and tritium /3 decay 
searches [61] as well as with the cosmological bound on the sum of light neutrino masses [62]. 
However, a positive signal of neutrinoless double beta decay in the next years as well as 
a cosmological measurement of the sum of neutrino masses of the order of 0.1 eV would 
certainly rule out this type of model. Therefore, future experimental improvements in the 
neutrino physics will be decisive for testing the viability of the SU(5) x Z4 model. 
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A Matter and Higgs representations 
Fermionic representations 

The fermionic fields in the model decompose in terms of the SM gauge quantum numbers 
as: 

5* = (1" e L, 

10 = Q e n'^ e e^ (A.i) 

24 = P8 © P3 © P(3,2) © P(3,2) © PO ■ 

The fermionic representations 5* and 10 can be written as 

5*=(cZ%, 5*=£,/, (A.2) 

and 

where a, P, j = 1, 2, 3 and i,j = 4, 5. The fermionic field p3, triplet of SU(2), belonging to 
the adjoint representation can be written as, 

1 / Ps' V2p,+\ 

where 

P3^ = '^^ Ps' = pI (A-5) 

Higgs representations 

The Higgs content of the model decomposes as 
5h = Ti® Hi, 

2Ah = Sg © S3 © S(3^2) © S(3_2) © ^0 , (A.6) 

A5h = 5(8,2) 3 © %,1) 1 ® -5(3,3) 1 © %,2) 7 ® %,1)4 ®T2® H2, 
in ~5 5 ~ID 5 
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where we have included for completeness the hypercharged properly normalised. The Hi 
and H2 are the usual Higgs doublets and Ti and T2 are the colour triplets. The 45 Higgs 
representation, which the explicit decomposition is given in ref. [27], obeys to the following 
relations, 

5 

45^ = -45i' and 1^457=0, (A.7) 

The different contributions for the beta coefficients bi of each extra particle besides the 
2HDM content are given in table 4. 



Table 4. Summary of the bi constants for relevant particles in the model. 



2HDM p3 ps jO(3,2) ^1,2 S3 Sg S(8,2) 

h 21/5 I 3^ I 

62- 3 |oi o|o I 

63- 7 02| |01 2 



B The Potential 

In this section we give explicitly the terms of the Higgs potential. Notice that index H on 
the Higgs fields is dropped in the following expressions. The potential V is divided into six 
parts as follows: 

V (5, 24, 45) =Vi (5) + V2 (24) + V3 (45) + V4 (24, 45) 

(B.l) 

+ 1/5(5,24) + ^6(5,45) , 
where each parcell are given by: 

V,{5) = -^5^5l + ^{5"5*J, (B.2a) 



V2 (24) = - ^ 24^ 24^ + ^ (24^ 24^ ) + ^ 24^ 24^ 24^ 

^ 2 V y d 2b) 

+ ^24^4^ 24] 24^, 
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^^3 (45) = - ^ 45f 45*^^ + A4 (45^45*^^) 

+ A545f 45*^^45f 45*;^^ + A645f 45*^^45^ 45*^^ 

+ A745f 45*^^45745*^5 + A845^^45*^,45S'^45*;^ (^-2^) 

+ A945^^45*?,45-45*i^ + Aio45^^45*?,45f 45*^, 

+ An45r45*^,45^^45*L, 
^4(24,45) =a2 45f 24^45*^^ + A12 45f 45*;^^ 24^24^ 

+ Ai3 45f 24^24^45*]^ + Au 45f 24;524f 45*^^ 

(B.2d) 

+ A15 45^24724^45*^5 + Ai645f 24^24^45*5;^ 
+ Ai7 45f 24:^24^45*^^, 
^5 (5, 24) = as 5^24^5^ + Aig 5^5^24^24^ + A19 5^24^24^5^ , (B.2e) 

(5, 45) = A20 45f 45*^.5|5'5 + A21 45f 5;45*;^.5^ + A22 45^45*^55^5^ . (B.2f) 
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